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The XMM-Newton view of NGC 6251 
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Abstract. We report on the nuclear X-ray properties of the radio galaxy NGC 6251 observed with XMM-Newton. 
NGC 6251 is a well-known radio galaxy with intermediate FRI/II radio properties. It is optically classified as a 
Seyfert 2 and hosts a supermassive black hole with mass ~ 6 x 10* Mq. The 0.4-10 keV EPIC pn continuum is 
best fitted by two thermal components {kT ~ 0.5 and 1.4 keV, respectively), plus a power law with photon index 
r ~ 1.9 absorbed by a column density A^h ~ 5 x 10^" cm~^. We confirm the previous ASCA detection of a strong 
iron line. The line, resolved in the EPIC pn spectrum, is adequately fitted with a broad (cr ~ 0.6 keV) Gaussian 
at rest-frame energy 6.4 keV with EW ~ 220 eV. We also detect, for the first time, short-term, low-amplitude 
variability of the nuclear flux on a timescale of a few ks. The spectral properties argue in favor of the presence 
of a standard accretion disk, ruling out the base of the jet as the sole origin of the X-rays. The moderate X-ray 
luminosity and lack of strong intrinsic absorption suggest that NGC 6251 is a "pure" type 2 AGN which lacks a 
broad-line region. 
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1. Introduction 

NGC 6251 (z=0.024) is a giant elliptical galaxy hosting a 
supermassive black hole with mass A/bh ~ 4 — 8 x 1O*M0 
(Ferrarese & Ford 1999), as measured with HST. At radio 
wavelengths, the source has been well studied. It shows 
the typical morphology of a FanarofF-Riley II (FR II) 
with a very narrow and long (> 4') jet extending N-W 
(e.g., Saunders et al. 1981), which emits at X-ray wave- 
lengths (Birkinshaw & Worrall 1993; Mack et al. 1997). 
However, based on its radio power at 178 MHz, NGC 6251 
is classified as an FR I (e.g., Owen & Laing 1989). Thus, 
NGC 6251 occupies an interesting niche in the study of 
radio galaxies as an intermediate object (from the radio 
perspective) between low- and high-power sources. 

In the optical, NGC 6251 has been classified as a 
Seyfert 2, based on the presence of permitted and for- 
bidden emission lines with FWHM ~ 600 km s~^ (e.g., 
Shuder & Osterbrock 1981). Previous X-ray studies of this 
galaxy with i^OS'^T' showed the presence of an unresolved 
nuclear source embedded in diffuse thermal emission as- 
sociated with the galaxy's halo (Birkinshaw & Worrall 
1993). ASCA observations showed that the nuclear spec- 
trum could be described as a moderately absorbed power 
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law continuum and also suggested the presence of a strong 
(EW'~ 400 eV), narrow line at ~6.7 keV (Sambruna et al. 
1999; Turner et al. 1999). 

Despite the intensive study of this source at all wave- 
lengths, the nature of the accretion process in NGC 6251 
is still a matter of debate. Based on the radio-to- X-ray 
spectral energy distribution, Ho (1999) suggested that an 
Advection-Dominated Accretion Flow (ADAF) is present 
in the nucleus of NGC 6251. On the other hand, Ferrarese 
& Ford (1999) and Melia et al. (2002) favored a standard 
accretion disk. Similarly, the origin of the X-ray emission 
from the nucleus is uncertain: Turner et al. (1997) ad- 
vocated a typical obscured Seyfert 2 spectrum based on 
the ASCA data, while Hardcastle & Worrall (2002) ar- 
gued in favor of an origin from the base of the relativis- 
tic jet for the soft X-rays. Further support for the jet- 
dominated hypothesis was recently claimed by Chiaberge 
et al. (2003) based on spectral energy distribution argu- 
ments. On the other hand, Guainazzi et al. (2003) favor a 
scenario with two main spectral components: a blazar-like 
spectrum dominating the high-flux state and a Seyfert-like 
spectrum emerging during the low-fiux state. 

With its larger sensitivity in the 0.3-10 keV band, the 
EPIC camera on-board XMM-Newton is an ideal instru- 
ment for investigating the origin of the X-rays and the 
nature of the accretion process in NGC 6251 through tim- 
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ing and spectral analyses. Motivated by the above con- 
siderations, we observed NGC 6251 with XMM-Newton 
in AOl for 50 ks. Here we report the results from the 
analysis of the nuclear light curve and spectrum. X-ray 
emission from the extended jet and halo will be dis- 
cussed elsewhere (Sambruna et al. 2003). Throughout 
the paper we use a Friedman cosmology with iJp = 
75 km s~^ Mpc~^ and go = 0.5. 

2. Observations and data reduction 

We observed NGC 6251 with XMM-Newton on 2002 
March 26 for a duration of ^ 41 ks with the EPIC pn, 
and for - 49 ks with EPIC MOSl and M0S2. All of the 
EPIC cameras were operated in full-frame mode with a 
medium filter for the MOS cameras and a thin filter for the 
pn. The recorded events were screened to remove known 
hot pixels and other data flagged as bad; only data with 
FLAG=0 were used. The data were processed using the lat- 
est CCD gain values, and only events corresponding to 
pattern 0-12 (singles, doubles, triples, and quadruples) in 
the MOS cameras and 0-4 (singles and doubles only, since 
the pn pixels are larger) in the pn camera were accepted. 
Arf and rmf files were created with the latest available re- 
lease of the XMM-Newton Science Analysis Software (SAS 
5.4). Investigation of the full-field light curves revealed 
the presence of a period of background flaring at the end 
of the observation. These events were excluded, reducing 
the effective total exposure time to ~ 36 ks for the EPIC 
pn and ~ 43 ks for the MOS cameras. The RGS data of 
NGC 6251 have a signal-to-noise ratio [S/N) that is too 
low for a meaningful analysis. Background spectra and 
light curves were extracted from source-free regions on the 
same chip as the source. There are no signs of pile-up in the 
pn or MOS cameras according to the SAS task epatplot. 
With an extraction radius of 35" the detected count rates 
in the energy range 0.4-10 keV are (0.519 ± 0.003) s^^ 
for the MOSl, (0.533 ± 0.004) s'^ for the M0S2, and 
(1.810 ±0.007) s~^ for the pn. For comparison with previ- 
ous broad-band X-ray satellites the count rates measured 
by the SISO on-board ASCA and the MFCS on-board 
BeppoSAX were ~ 0.07 s~^ and ^ 0.06 s~^, respectively 
(Sambruna et al. 1999; Guainazzi et al. 2003). 

Inspection of an archival Chandra observation of 
NGC 6251 indicates that no serendipitous sources are 
present within the EPIC extraction radius and that 
the contribution from the resolved kpc jet is negligible. 
However, no entirely reliable information on the nuclear 
properties of NGC 6251 can be drawn from Chandra data, 
due to the unfortunate location of the nucleus on the 
CCDs' gaps and to the pile-up. 

The EPIC spectra were rebinned such that each spec- 
tral bin contains at least 100 counts for the pn and 40 
counts for the MOS cameras in order to apply mini- 
mization, and fitted jointly using the XSPEC v . 11 . 2 soft- 
ware package (Arnaud 1996). The quoted errors on the 
derived best-fitting model parameters correspond to a 



2.1 5 r 




4 



Time (s) 

Fig. 1. EPIC pn+M0Sl+M0S2 light curves of the 
background-subtracted count rate in the 0.8-10 keV band 
(top panel) and of the X-ray color 1.8-10 keV/0.8-1.5 keV 
(bottom panel). The extraction radius is 35"; time bins are 
3000 s. The dashed lines indicate the average values. 

90% confidence level for one parameter of interest (i.e., 
a Ax^ — 2.7 criterion) unless otherwise stated. 

3. EPIC light curve of the nucleus 

We first study the 0.3-10 keV light curve obtained with 
the EPIC pn, which is the most sensitive instrument on- 
board XMM-Newton, using an extraction radius of 35"and 
a time bin of 3000 s (changing the extraction radius or the 
time bin has a negligible impact on our results). At ener- 
gies below 0.8 keV no statistically significant variability is 
detected according to a x^ test. This might be partly due 
to a combination of dilution from the thermal emission 
and intrinsic absorption. However, the small intrinsic ab- 
sorption and the relatively small contribution of the ther- 
mal component in the low energy range (see §4) suggest 
that the very-soft (0.4-0.7 keV) flux might be intrinsi- 
cally constant or, at most, variable at a very low level. We 
therefore focus on the 0.8-10 keV energy band. The com- 
bination of pn, MOSl, and M0S2 light curves increases 
the statistical significance of the variability. Figure ^ (top 
panel) shows the 0.8-10 keV background-subtracted EPIC 
pn-|-M0Sl-|-M0S2 light curve of NGC 6251. The bottom 
panel shows a plot of the hardness ratio, defined as the ra- 
tio of the 1.8-10 keV count rate to the 0.8-1.5 keV count 
rate, versus time. Low-amplitude flux variability is present 
on timescales of a few ks. A test yields a probability 
that the count rate is constant of P^2=0.9%. The back- 
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ground count rate is as low as 5% of the average source 
count rate. The variability of the hardness ratio is not 
significant (P^2 -22%). 



Table 1. Short-term X-ray variability of NGC 6251 



Energy band 


(r) 


Xrcd 




P a 
var 


(keV) 




(14 d.o.f.) 


(%) 


(%) 


0.4-0.7 keV 


0.49 


0.54 


89.6 




0.8-1 keV 


0.34 


1.40 


14.8 


2.2 ± 1.4 


1.2-1.5 keV 


0.36 


1.48 


11.7 


2.1 ± 1.3 


1.8-3 keV 


0.41 


1.59 


8.1 


2.3 ± 1.2 


3-10 keV 


0.39 


1.66 


6.2 


2.8 ± 1.3 



The errors on Fvar are calculated as in Edelson et al. (2002) 
and should be considered conservative estimates of the true 
uncertainty. 



The study of energy-dependent variability is hampered 
by the presence of photons emitted from the diffuse com- 
ponent that dilutes the intrinsic variability of the source, 
especially at low energies. Therefore, we have 1) divided 
the 0.8-10 keV range into several sub-bands; 2) checked 
with a test the presence of variability in each sub-band; 
3) merged contiguous variable sub-bands, creating four fi- 
nal bands with roughly the same count rates. The energy- 
dependent variability has been characterized by means of 
the chance probability and the fractional variability pa- 
rameter, Fvar — (cr^ — A^)^/^/ (r), where is the variance, 
(r) the unweighted mean count rate, and the mean 
square value of the uncertainties associated with each in- 
dividual count rate. The results are reported in Tabled 
For comparison, the 0.4-0.7 keV results are also reported. 
No clear variability trend with energy is found. 

NGC 6251 has been observed by several X-ray satel- 
lites over the last decade. Therefore, its long-term tempo- 
ral behavior can be investigated. A comparison of the 2- 
10 keV absorbed flux detected by ASCA in October 1994 
(Sambruna et al. 1999), Fasca ^ 1.4 x 10' erg cm'^ 
s""'^, with values more recently measured with BeppoSAX 
(observation carried out in July 2001), Fsax ~ 4.7 x 
10~^^ erg cm~^ s~^, and with XMM-Newton (see below) 
^XMM ~ 4 X 10"^^ erg cm~^ s~^, suggests an increase by a 
factor of three over the last nine years. Although this con- 
clusion might be weakened by systematic errors related to 
the instrumental cross-calibration uncertainties, the spec- 
tral uncertainties, and the different extraction areas, the 
substantial flux difference measured argues in favor of gen- 
uine variability. Observations with the same satellite are 
needed to carefully investigate the long-term variability. 
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Fig. 2. EPIC pn unfolded spectrum of NGC 6251 de- 
scribed by a thermal component plus an absorbed power 
law. The lower panels show the data-to-model ratio for the 
pn, MOSl, and M0S2 cameras, respectively. The edge at 
0.53 keV in the top panel is due to neutral oxygen in the 
absorbing material. 



4. The nuclear spectrum 

4.1. The continuum 

Motivated by the controversial results of previous X-ray 
spectral studies, we fitted the EPIC spectra in the 0.4- 
10 keV range, where the instruments are best calibrated, 
with models of increasing complexity. To account for resid- 
ual calibration problems between the EPIC cameras, we 
first fitted separately the pn, MOSl, and M0S2 spec- 
tra. All the spectral models assume a column density, 
obscuring all components, fixed at the Galactic value of 
Nn = 5.65 X 10^° cm'^ (Dickey & Lockman 1990). In some 
models an additional intrinsic absorber at the redshift 
of NGC 6251 is included. A power-law model modified 
by Galactic absorption gives unacceptable fits: /d.o.f . 
are 927.8/443, 401.2/232, 384.5/233, for the pn, MOSl, 
and M0S2, respectively. A first significant improvement 
(Ax^ ^ 80 for three additional parameters) is obtained 
by adding a soft component modeled as emission by a 
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Table 2. Spectral Model Parameters for NGC 6251. Parameters without errors have been fixed at the best-fit values 



a) EPIC continuum parameters 







thermal 






absorbed PL 








Camera 


kT 
keV 


Z 

Zq 


normapoc 
10-*cm-'5 


lO^Ocm-^ 


r 


normpow 

10~^ J'lkcV 




Count rate 


MOSl 




0.19 




0-0-0.9 


l-00_(, 04 


1 in+o.05 

'-■'-''-0.04 


243.8/229 


0.519 ± 0.003 


M0S2 


0.52l»;- 


0.18 




7 n+0.7 


^•y^-0.05 


1 ni+0.03 

^•■^^-0.06 


236.0/230 


0.533 ± 0.004 


FN 




0.18 




5.01°:? 


1 01+003 


1 OO+0.02 
J^-50-0.05 


548.4/440 


1.810 ±0.007 



thermali 



b) Best-fitting continuum: EPIC pn 

thcrmal2 



absorbed PL 



feTi 
keV 



Zi 
Zq 



keV 



Z2 
Zq 



norm 
10 



apec2 

*cm-5 



102°cm- 



10-' FiV 



-1 



1.36 



+0.34 



0.31 



2.25 



+ 1.37 



0.541 



0.17 



3.23 



+0.72 



5.1; 



-1.1 

-0.9 



1.91 



+0.08 
0.05 



1.91 



+0.08 
0.05 



520.3/435 



c) Line parameters: EPIC pn 



E 
keV 


a 
keV 


flux 
lO"'^ s"^ cm-2 


EW 
eV 


6 42+0-52 

"-^^-0.41 


0.58tr2i 


S c;+8-3 
8-0-3.8 


990+219 

^^■^-99 



* photon keV ^ cm ^ s 



collisionally-ionized plasma (apec in XSPEC; Smith et al. 
2001). However, a fair representation of the EPIC spec- 
tra is only obtained with a thermal component plus an 
absorbed power law. The column density of cold gas to- 
ward the nucleus TVh ~ 5 x lO'^" cm"^, although small, 
is statistically highly significant: Ax^ ~ 70 — 100 for one 
additional parameter. We have investigated the possibility 
that the intrinsic absorber is ionized (using the absori in 
XSPEC) or partially covering the nuclear source (using the 
zpcf abs model in XSPEC). However, neither of these mod- 
els resulted in a significant improvement of the fit. The 
best-fit parameters for the model consisting of a thermal 
component plus an absorbed power law, are reported in 
Table 2a (the abundances were frozen at the best-fit val- 
ues) . Figure [S] shows the unfolded spectrum of the EPIC 
pn and the data-to-model ratio for the three EPIC cam- 
eras. A close inspection of residuals shows a more complex 
pattern in the pn data. Indeed the pn fit is further im- 
proved by adding a second thermal component. This hot- 
ter thermal component is consistent with the one detected 
in the circumnuclear gas around NGC 6251 (Sambruna et 
al. 2003). On the other hand, the fits to the MOS spec- 
tra are not improved significantly because of their lower 
S/N . We also tried an alternative model for the pn con- 
tinuum with two power laws plus a thermal component. 
However, the resulting spectral fit, {kT = 0.59io;o5 keV, 
Z ^ 0.18 Zq, Pi = 2.68115^3, ^2 = l-GStaos) is sta- 
tistically worse than the one obtained with two thermal 



models and one absorbed power law (A^^ = -1-9.6 for the 
same d.o.f.). 

In addition, the residuals of the pn and M0S2 data 
show a clear excess in the 6-7 keV range, suggesting the 
presence of a line (see §4.2). The best- fitting continuum 
parameters are summarized in Table 2b and the unfolded 
spectrum plus data-to- model ratio is shown in Figure |2| 

The total absorbed fluxes from 0.4-2 and 2-10 keV are 
~ 2.6 X 10^^^ erg cm^^ s^^ and ^ 4.0 x 10^^^ erg cm^^ 
s"^, respectively, corresponding to intrinsic (absorption- 
corrected) luminosities of Lo.4_2 koV ~ 4.4 x 10"*^ erg 
s^^ and L2_io kcV ^ 4.8 x lO"*^ erg . The contri- 
bution of the thermal component is ~ 10% from 0.4- 
2 keV and negligible above 2 keV. This leads to an X-ray 
luminosity associated with the power-law component of 
-^0 4-10 keV ^ ^-^ ^ ^^^^ "^^S s""^. This value can be further 
increased up to 1.7 x lO'*'^ erg s"^ by extending the en- 
ergy range up to 100 keV. In fact, analysis of a BeppoSAX 
observation of NGC 6251 shows that the PDS detects the 
source up to energies higher than 100 keV (Guainazzi et 
al. 2003). 

4.2. The Fe line 

To investigate the possible presence of an Fe Ka line, we 
used only pn data. In the MOS cameras, the rapid decrease 
of quantum efficiency with increasing energy reduces the 
S/N above 6 keV, hampering a thorough analysis of the 
Fe Ka range. To reduce the influence of the complex con- 
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Fig. 3. EPIC pn unfolded spectrum of NGC 6251 
and data-to-model ratio. The spectral model is 
wabs(apec-|-apec-|-zwabs(powerlaw-|-zgauss)). 




Gaussian line included 



2 5 1C 

E (keV) 

Fig. 4. Top panel: EPIC pn data-to-model ratio from 2~ 
10 keV of the best-fit continuum model for the EPIC pn 
data. Bottom panel: same as above but with a Gaussian 
line included. Energies are in the observer's frame. 

tinuum modeling on the line detection and properties, we 
restrict the energy range to 2-10 keV. The 5-7.5 keV band 
was initially excluded to allow a better determination of 
the photon index. A simple power law model produces a 
satisfactory fit to the the 2-5 keV and 7.5-10 keV win- 
dows of the pn spectrum, with = 134.1 for 123 degrees 
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of freedom. Including the 5-7.5 keV energy range, the fit is 
still formally acceptable (x^ — 173.1, 158 d.o.f). However, 
residuals around 6-7 keV (Fig. 0] top panel) suggest the 
presence of a broad line. Indeed, adding a Gaussian line 
improves the fit significantly (Ax^=14 for 3 additional 
parameters), as can be seen from the residuals in the bot- 
tom panel of Fig2| For testing purposes, we applied the 
same procedure to M0S2 data. A marginally significant 
(Ax^ = 7.1 for three additional parameters) line at ~ 6.4 
keV is detected, with parameters consistent within the er- 
rors with the ones inferred using pn data. Using the entire 
0.4-10 keV range and the best-fitting model, the statisti- 
cal significance of the line detection is further enhanced: 
Ax^ ^ 15.5 for both pn and M0S2. 

We also tried more complex fits to the broad line using 
specific disk models in XSPEC, such as diskline and laor. 
The former model parameterizes the expected line profile 
from a disk around a Schwarzschild black hole (Fabian et 
al. 1989), whereas the latter describes the line profile in 
the case of a rotating Kerr black hole (Laor 1991). In both 
cases, the fitted inner radius of the disk is located within 
a few Schwarzschild radii from the black hole; however, 
due to the limited S/N, the line parameters are poorly 
constrained. 

In summary, we find evidence for a broad line at rest- 
frame energy ^ 6.4 keV with EW~ 220 eV and width 
FWHM~ 64,000 km s~^, which we interpret as the fluo- 
rescent Fe Ka line commonly observed in Seyfert galaxies. 
The line parameters with the 90% confidence errors are 
given in Table 2c. 

5. Discussion and Conclusions 

A detailed discussion of the origin of the X-rays from weak 
radio galaxies has been reported by Gliozzi et al. (2003) 
for the specific case of NGC 4261. For that object, a jet 
origin, claimed on the basis of a tight correlation observed 
between radio and X-ray core luminosities, was disputed 
on the basis of variability properties and energetic con- 
siderations. On the other hand, for NGC 6251 rapid vari- 
ability is detected but no clear flux or spectral variability 
trend is found, and the investigation of the origin of the 
X-rays must rely upon the time-averaged spectral results. 

The strongest argument against a jet origin comes from 
the detection of a prominent broad Fe Ka line. The large 
EW is consistent with the average value found in Seyfert 1 
galaxies (e.g., Nandra et al. 1997) and with the upper lim- 
its found in brighter radio- loud AGN (e.g., Sambruna et 
al. 1999), and is much larger than typical values found 
in broad-line radio galaxies observed with RXTE (e.g., 
Eracleous et al. 2000), arguing against a significant con- 
tribution from the jet. A Doppler-broadened line with 
FWHM~ 64, 000 km s^-^ can only be produced in the in- 
ner region of an accretion flow by fluorescence (although a 
deeper X-ray exposure is necessary to constrain the inner 
radius of the line-emitting region). Contributions from the 
BLR are excluded here because the source lacks optical 
broad lines (e.g., Shuder & Osterbrock 1981). Alternative 
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scenarios invoking a jet with a wide opening angle illu- 
minating the putative obscuring torus can only produce 
narrow emission lines (Wozniak et al. 1998). 

In addition, the well constrained photon index, F ~ 
1.9, agrees very well with the average value found in 
Seyfert 1 galaxies (e.g., Nandra et al. 1997); however, it is 
not inconsistent with photon indices derived for X-ray jets 
of FR I radio galaxies, even though the latter are generally 
poorly constrained (e.g., Hardcastle et al. 2002). 

The small absorbing column is in good agreement 
with the low visual extinction Ay ~ 0.6 inferred for the 
dusty disk (Ferrarese & Ford 1999), assuming the Galactic 
gas-to-dust ratio. The lack of strong intrinsic absorption 
in NGC 6251 argues against obscuration by a canonical 
molecular torus (see, e.g., Chiaberge et al. 2002). The issue 
of intrinsic absorption in FR I galax;ies will be addressed 
in a forthcoming paper (Donato et al. 2003, in prep). 

Assuming that (most of) the X-rays are not associated 
with the jet, we investigate the nature of the accretion 
flow. The main question is whether an ADAF or a stan- 
dard disk model provides a better description of the accre- 
tion process in NGC 6251. To answer this question, we first 
note that the ADAF interpretation proposed by Ho (1999) 
was mainly prompted by the low value inferred for the 
ratio iboi/-^^Edd, where Lboi = 8.2 x 10''^ erg s~^ was ob- 
tained by integrating the spectral energy distribution from 
radio to X-rays. However, because of the paucity of data 
points, the estimate of Lboi is highly uncertain. Indeed, 
Ferrarese & Ford (1999) and Meha et al. (2002) indepen- 
dently derive a value of Lboi larger by a factor of a few with 
respect to the value derived by Ho. Direct evidence that 
the value derived by Ho underestimates Lboi comes from 
the value of X-ray luminosity Lo.4-10 keV ~ ^-^ ^ ^^^^ 

associated with the power-law component. Therefore, 
assuming that Lq\_-^q ^.^y represents '~30% of Lboi (see, 
e.g., Elvis et al. 1994), we obtain a more reasonable but 
still conservative estimate of Lboi of the order of 3 x 10''^ 
erg s^^, consistent with the estimate given by Ferrarese 
& Ford (1999). However, the modest value of the ratio 
Lboi/'t'Edd derived from this assumption does not allow us 
to discriminate between the competing accretion scenar- 
ios. 

An independent argument favoring a standard accre- 
tion disk is the detection of a broad Fe Ka line. As already 
mentioned, the most natural explanation for the broad line 
detected at ~ 6.4 keV is that it is a fluorescent Fe Ka line 
produced within the inner gravitational radii of a optically 
thick accretion disk. 

Direct evidence against radiatively ineflicient solutions 
comes from the estimate of the radiative efliciency rj, 
which can be derived comparing the accretion luminosity 
-t-accr = V^siccrC^ to Lboi- A rough estimate of the accre- 
tion rate is given by the Bondi value, Meondi = '^t^RaPaCs 
where i?A — GM/c^ is the accretion radius, Cg the sound 
speed, and pA the density at the accretion radius. Since 
the soft thermal component (kT ~ 0.55 keV) is spa- 
tially unresolved with XMM-Newton (Sambruna et al. 
2003), it is appropriate to choose the instrument reso- 



lution, 10", as the angular size of the thermal compo- 
nent. With this assumption, and considering normal cos- 
mic composition {ue/up = 1.18), the emission measure 
given by the normalization of the soft thermal component, 
[10-iV(47r(DA(l + z))^)] J UeUHdV ~ 3 X lO"'' cm'^ 
(where Da is the angular size distance to the source in 
cm, and rie and njj are the electron and H densities in 
cm~^) yields a proton density of Up ~ 0.06 cm~^, which 
translates into Msondi — 8 x 10"'^ Mq yr^^. A direct com- 
parison of this value with the estimated Lbd leads to a 
radiative efficiency of the order of 77 ^ 0.06, fully con- 
sistent with the standard disk scenario. This conclusion 
is supportc^d by tlic^ striking difference^ observed between 
NGC 6251, that radiates at the Bondi limit, and low lu- 
minosity AGN observed with Chandra (Loewenstein et al. 
2001) which show Lboi <C Lsondi and are interpreted in the 
ADAF framework. 

A final important issue to be addressed is the nature 
of the Seyfert nucleus in NGC 6251. Even though histori- 
cally this source has been classified as a type 2 AGN, its 
X-ray properties as seen with XMM-Newton seem to be 
more typical for a type 1 AGN. Interestingly, we note that 
the value inferred for Lboi/LEdd translates into M / M-^dd 
in the range where recent models (e.g., Nicastro 2000, 
Laor 2003) predict that a broad-line region (BLR) can- 
not exist. The absence of broad optical lines and the low 
value inferred for the cold absorption from the X-ray spec- 
tral analysis support the hypothesis that NGC 6251 is 
a "pure" type 2 AGN, i.e., without BLR. This resuh is 
consistent with the findings of Steffen et al. (2003), who, 
studying sources of the X-ray background, find that the 
fraction of objects with broad lines drops sharply below 
Lx ~ 10^'' erg s~^. An alternative interpretation might be 
that the source is Compton-thick and the observed X-rays 
are only the scattered component. A possible way to test 
this hypothesis is based on the location of NGC 6251 in the 
Fe EW - Lx/Loiii plane (Bassani et al. 1999). However, 
such analysis does not provide any firm conclusion, since 
NGC 6251 is located half way between the Compton thick 
region and the Seyfert 1 region. The rapid X-ray variabil- 
ity, the rather weak iron line for a type 2 AGN, and the 
BeppoSAX results (Guainazzi et al. 2003) seem to argue 
against the Compton-thick hypothesis. However, deep ex- 
posures at hard X-rays are necessary to discriminate be- 
tween the two competing scenarios. 

A recent cross-correlation study of the FIRST and 2DF 
catalogs (Magliocchetti ct al. 2002) finds that the ma- 
jority of radio sources do not show any optical emission 
lines, while a minority have spectra similar to NGC 6251. 
If NGC 6251 with its Seyfert 1-like nuclear properties is 
representative of the whole group, it follows that there is 
a significant fraction of sources which do not fit the stan- 
dard view of AGN. Detailed multiwavelength studies of 
their nuclear properties are needed. 
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